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Abstract 

In the framework of a two Higgs doublet model type I and type II, we calculate limits on neutrino masses for the 
different types of neutrinos, by using the experimental bounds on their magnetic dipole moments. This is carried out 
by analyzing diagrams of Cherenkov neutrino decays with a charged Higgs into the loop, coming from the two Higgs 
doublet model (2HDM). Such constraints are translated into allowed regions in the free parameters of the models, for 
each neutrino flavor. 

The analysis was performed by sweeping the charged Higgs mass between (100 — 9Q0)GeV and taking into account 
the experimental constraints for tan in the 2HDM type I and II, obtaining contributions close to the experimental 
thresholds for muon and tau neutrinos, while for electron neutrino the relevant contribution comes from standard 
model and keeps out of the reach of forthcoming experiments. 

PACS: 41.20.Cv, 02.10.Yn, 01.40.Fk, 01.40.gb, 02.30.Tb 


1 Introduction 

The phenomenon of neutrino oscillations, has been supported by the discovery of flavor conversions of neutrinos from 
different sources, like the atmospheric neutrinos made by Super-Kamiokande in 1998 [J, or more recently by T2K 
Collaboration]^. These oscillations occur among at least three types of flavors of neutrinos: electron, muon and tau 
neutrinos. 

K) = E U °k I"*) ( Q = e ’ /*> T ) (*) 

k 

all neutrinos produced and observed so far, have left-handed helicities, while all antineutrinos have right-handed helicities. 
Neutrino oscillations provided the first glimpse of physics beyond the standard model of particles. Now, since neutrino 
oscillations are sensitive only to the difference in the squares of their masses, such a phenomenon requires that at least 
two neutrino species have nonzero mass. The transition probability between different flavors can be approximated by 

—, ( Am? -L\ 

M = E U: k U 0 k U aj U - Bj exp (2) 

where Pv a -n, g ( t ) is the probability that after travelling a distance L , a neutrino with flavor i/ Q converts into a neutrino 
with flavor i/p. As for the neutrino masses, we have only upper boimds hitherto[3] 

m Ue < 2.0beV , m„„ < 0.19MeU , m Vr < 18.2Mel/ (3) 

All elementary fermions in the Standard Model are Dirac fermions. Nevertheless, the nature of the neutrino is not 
yet definitely settled and depending on the model the neutrino can be either a Majorana or Dirac fermion. On the 
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other hand, despite the neutrinos do not carry electric charge, they can participate in electromagnetic interactions by 
coupling with photons via loop diagrams, and like other particles the electromagnetic properties can be described by 
electromagnetic form factors (EFF’s). For example, by means of its multipole moments, neutrinos can be sensitive to 
intense electromagnetic fields, and such intense fields can exist in nature. It has been suggested that there could be 
sources of magnetic fields of order (10 13 — 10 18 ) G , as it could be the case during a supernova explosion or in the vicinity 
of special groups of neutron stars known as magnetars|4]. 

On the other hand, present limits on the scalar sector in the standard model, still permits the possibility of an extended 
Higgs sector. We shall study one of the simplest extension of the scalar sector of the standard model, the so-called Two 
Higgs Doublet Model (2HDM) in which we add a second Higgs doublet with the same quant urn numbers of the first. 
There are many motivations for this model, one of this is the fact that the SM is unable to generate a baryon asymmetry 
of the universe of sufficient size, or to explain the mass hierarchy in the third generation of quarks. Two Higgs Doublet 
models are possible scenarios to solve these problems, due to the flexibility of their scalar mass spectrum and the existence 
of additional sources of CP violation. In addition in the Minimal Supersymmetric Standard Model (MSSM), a second 
doublet should be added in order to cancel anomalies [5]. 

The coupling of neutrinos with photons occur via loop diagrams. In the Standard Model (SM), the loop corrections 
have the form of vertex diagrams and vacuum polarization diagrams. When a second doublet of scalars is included in the 
spectrum, further corrections appear by replacing the vector bosons W ± by charged Higgs bosons H ± . Our goal is to 
characterize the corrections to the EFF’s coming from the new physics, and particularly on the region of parameters in 
which such factors become near the threshold of detection. In the region of parameters in which the threshold of detection 
is reached, we obtain bounds on neutrino masses. 

The structure of this paper is as follows. In section 0 we discuss briefly the implementation of neutrino Dirac masses 
in SM. In section (3] we discuss the general form of the EFF’s for neutrinos. In section U we describe briefly the two Higgs 
doublet Model (2HDM), particularly the models of type I and of type II as well as the implementation of neutrinos masses 
into those models. In section 0 we characterize the loop diagrams coming from the 2HDM that contributes to the EFF’s 
of the neutrino. In section 0 we find upper bounds for the neutrino masses in the framework of the 2HDM type I and 
II, by using the allowed values of the free parameters of the model, as well as the experimental limits for the magnetic 
dipole moments of such neutrinos. Finally, section 0 yields our conclusions. 


2 Neutrino Dirac mass term in SM 


There are several ways to incorporate neutrino masses within the SM or its extensions, in order to explain the observed 
neutrino oscillations. We shall use a simple form which consists of adding right-handed singlets of neutrinos fields (&V/0 
corresponding to each charged lepton. This insertion implies new gauge invariant interactions in the Yukawa sector 

3 

- ^Yukawa = Y Y fo^aL^a'R + h.C (4) 

Q=C,/i,T a'=l 


where / a ,a' is a matrix with new coupling constants, Tp a L is the left-handed lepton doublet and $ is the SM Higgs doublet, 
with 4* = 



(5) 


a nonzero vacuum expectation value (VEV) of the Higgs doublet induces the spontaneous symmetry breaking from 
SU (2) l x U (l)y to U (1)q- In turn, the VEV also provides the neutrino Dirac mass term 




_ y 

y/2 ^ 

v Q=e,n, 


Y {a^cL^c'R + h.C 


a =1 


( 6 ) 


In general the matrix / Q , Q ' is a complex 3x3 matrix, the massive neutrino fields are obtained through the diagonalization 
of CD maaa > this can be done by diagonalizing / a ,a' with a biunitary transformation 

(P* / V) k j = m k S kj 


(7) 
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where m* is a diagonal matrix with real and positive values. 

The flavor eigenstates z/ c , and v T are superpositions of the mass eigenstates 

3 3 

v aL = ^2 ^o,k^kL ; Vq'R = ^2 Va',j u jR 

A:=l j= 1 

or equivalently, the mass eigenstates i/^ are mixtures of flavor eigenstates. The mixing matrix U is called the PMNS 
matrix due to Pontecorvo|B], Maki, Nakagawa and Sakata[7]. 

The result of the diagonalization gives Dirac mass terms of the form 

3 

- C D ma ,. = Y, ' n ^k ( 8 ) 

A:=l 

with the Dirac fields of massive neutrinos given by 

Vk = I'kL + VkR (9) 


3 The electromagnetic form factors (EFF’s) 



Figure 1: Effective coupling of two neutrinos with a photon 

(u (p, A) \J EM (*) \u (p', A')) = 77 (p, A) (l, q) u (p 7 , A') (10) 

To find all the EFF’s, we use the general expression for the current @,0, where q^ = p' ll — p,i, = pj, + Pn are the 
four-momenta shown in Fig. [T] and u (p, A), u (p 7 , A 7 ) are the initial and final fermion states respectively. Further, A p are 
matrices of couplings acting on the spinors. The matrices A^ have some interesting properties 

• The first condition is that the arrangement A^ must be a 4-vector, i.e. must be Lorentz covariant. 

• The second condition is hermiticity of the associated current, i.e., J^ EM = J EM which implies 

A„ ((,</) = 7°A»(i, -9)7° (11) 

• The current conservation or gauge invariance f = 0 can be recast into 

q^u (p 7 , A 7 ) A m (l,q) u (p, A) = 0 (12) 

The most general expression for A M (l, q) reads [ID] 

A /1 (?) = f Q (q 2 ) 7 m + [ F M (q 2 ) i + Fe (q 2 ) 75 ] Vfivq 1 ' + Fa (q 2 ) (q 2 7 m ~ 9 m ^) 75 ( 13 ) 


where Fq , Fm, Fe and Fa represent the electric charge, dipole magnetic moment, dipole electric moment and anapole 
moment respectively. 
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The EFF’s show us how the particles are coupled with the photon at the tree level or in loop corrections. At the tree 
level we got the electric charge and one part of the contribution coining from the magnetic dipole moment. Now, if we 
consider the interaction with an external field A^ xt in the form 

C' X t = -eA* xt J™ (14) 

the so-called anomalous magnetic moment arises. Even uncharged particles may have magnetic dipole moment. However, 
for uncharged particles all dipole moments only appear in loop corrections. Just like the anomalous magnetic moment, 
the dipole electric moment and the anapole moment can be non-zero even for an uncharged particle[TT], We summarize 
some electromagnetic properties of charged leptons in table E 


l 

Mass(A''--v / ) 

MDM 

EDM (^) 

e 

0.51 

1.159 x 10 *hb 

Clx 10~ r ° 

t 1 

105.658 

1.159 x 

< 2 x 10"* 

T 

1.776 x 10° 

1.159 x 10 

2m T 

< 2 x 10" z 


dipole moment and EDM electric 

dipole moment 

Like other particles, neutrinos can be described by EFF’s with vertex functions. For neutrinos the magnetic and 
electric dipole moments are expected to be very small since they are likely proportional to the neutrino masses. For the 
anomalous magnetic moment the leading contribution is [IT] 


3 GFm Ui 

a Ui =- _ _ m. 




(15) 


Consequently, the neutrino magnetic moment is [15] 


H Ui = -—dviSvi 


(16) 




where /ib is the Bohr’s magneton. If neutrino couples to photons via such moments, the neutrino electromagnetic 
properties can be used to distinguish Majorana and Dirac neutrinos. For Dirac neutrinos the most relevant moment is 
Fm, because the other terms vanish in a CP —conserving scenario with an hermitian J ^ M , and are highly supressed owing 
to the soft violation of CP. On the other hand, for Majorana neutrinos only Fa is possible, because the other terms 
vanish owing to the self-conjugate nature of Majorana neutrinos. Table [2] summarizes the electromagnetic properties of 
massive neutrinos [H] 


/ 

Massl*" 1 ') 

Magnetic dipole moment 

i/ e 

< 2.2 x 10" b 

< 10.8 x 10 


< 0.17 

< 7.4 x 10“ 1U // B 

V r 

< 16 

< 5.4 x 10 ~‘fiB 


Table 2: Electromagnetic properties of massive neutrinos. 


4 The two Higgs doublet model with massive neutrinos 

The 2HDM contains five Higgs bosons in its spectrum]!^. The symmetry breaking is implemented by introducing a new 
scalar doublet with the same quantum numbers of the first one[l3]. In a CP-conserving scenario, the Higgs sector consists 
of: Two Higgs CP-even scalars (//°,h°) , one CP-odd scalar (A 0 ) and two charged Higgs bosons (H ± ). A key parameter 
of the model is the ratio between the vacuum expectation values 

tan 0 = — 

i>i 


( 17 ) 
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where v\ and V 2 are the vacuum expectation values of the Higgs doublets[13]. with values of 0 < 0 < §. 

The most general gauge invariant Lagrangian that couples the Higgs fields to leptons (with massless neutrinos) reads 



where $ 1,2 represents the Higgs doublets, and # 1,2 = io 2 ^ 1 , 2 * The superscript “0” indicates that the fields are not mass 
eigenstates yet, rjij and fij are non diagonal 3x3 matrices with (i, j) denoting family indices. Ej R denotes the three 
charged leptons and f iL denotes the lepton weak isospin left-handed doublets. 

It is customary to implement a discrete symmetry in the 2HDM in order to suppress some processes such as the Flavor 
Changing neutral currents (FCNC). In particular by demanding the discrete symmetry 


$1 — y <J>! ; $2 -► —$2 

D jR -y =f Dj R ; Ujr -> -Ujr (18) 

such kind of processes are efiminated at the tree-level. Here D jR and Ujr denote right-handed singlets of the down and 
up types of fermions. 


4.1 The 2HDM type I 

By taking DjR —> —DjR we arrive to the so-called 2HDM of type I. In this scenario, only $2 couples in the Yukawa sector 
and gives masses to all fermions. The Lepton Yukawa Lagrangian becomes 

-Cy = vZ’Xl^r + + h.c. (19) 

and the term of charged current of the Lagrangian with leptons yields 

-C Y = - M^Pl) vH + + h.c. (20) 


4.2 The 2HDM type II 

If we use DjR —> DjR we obtain the so-called 2HDM of type H. In this model $1 couples and gives masses to the down 
sector, while <f> 2 couples and gives masses to the up sector. In consequence, the Lepton Yukawa Lagrangian (with massless 
neutrinos) becomes 

-C Y = r#°7° L 5iE° fl + + h.c. (21) 

and the term of charged current of the lagrangian with leptons gives 


Cy = 


s/2M 


w 


[ (cot QM^Pl + tan &Mf ia9 Pnj IH + + h.c. 


An interesting aspect is that the limits for the parameter space (m//+,tan/3), for model type II are very similar to 
those obtained by considering the minimal supersymmetric scenario. 


4.3 2HDM type I and II with massive neutrinos 

The term J2J) inserted in the Yukawa sector of the standard model should also be inserted in the two Higgs doublet model 
for each doublet 4>,-. Nevertheless, when we implement the discrete symmetry iTTsl in a Lagrangian of the form Ql we 
observe that the term involving the doublet cannot appear, and that the extra term is the same in either model type 
I or model type H. 
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v v 

Figure 2: Loop corrections with leptons and W ± vector bosons in SM. 

.v 



Figure 3: Vacuum polarization with W ± vector bosons, and fermions denoted by / in SM. 



v v 

Figure 4: Loop corrections with leptons and H ± in the 2HDM. 



Figure 5: Vacuum polarization with H ± in the 2HDM. 

5 Radiative corrections in 2HDM 


The diagrams that contribute to the neutrino electromagnetic vertex in SM are displayed in Fig. [2] 

And for the vacuum polarization they are shown in Fig. [3] Within the framework of a 2HDM with massive neutrinos, 
we should add three new types of diagrams: two vertex corrections illustrated in Fig. (U and one correction to the vacuum 
polarization displayed in Fig. \5\ They arise by replacing W ± by H ± in the SM diagrams. 

Fig. [Hf) shows the vertex correction involving two charged Higgs bosons and one charged lepton into the loop 
(2H ± 1L). For this diagram, the general form of the contribution can be written as 




±1 L 


(9.0 = -* / 


d 4 k 

(2^ 


(2 k° + pg + p?) (aP L + bP R ) (# + m t ) ( cP L 
[(fc + pi)' 


dP R ) 


m H± 


[(* + P 2) 2 ~ m 2 H± j (k 2 - mf) 


where a, b,c and d are constants associated with the Feynman rules of the 2HDM, with Pr,l = (l± 7 5 ) /2. The 
contribution to the EFF’s of this diagram is described in the appendix, and in particular, the contribution to the magnetic 
dipole moment (MDM) is given by 


1 X 

A!h*il( 9 . 0 mdm = 1^2 J d * J d y - p?2 [( m " 

0 0 



+ (nn — 3 m u ) x + 2 x 2 m l/ 


i0 QfA qn [(ac + bd) + (bd - ac) 75 ] 


On the other hand, the diagram in Fig. [lie) with two leptons and one charged Higgs into the loop (2L1H ± ), gives a 
contribution of the form 

f d*k \ {.aPL+bP R ){t + p i + mi M* + + m,) ( cP L + dP R ) 1 

J (2tt) 4 j + Pl ) 2 _ m 2 j ^ + p 2 ) 2 _ m 2 j (^.2 _ m 2 H± ) j 
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from which we obtain the contribution of this diagram to the MDM, that is given by 

1 X 

A 2LIH± (?> 1 )mdm = J dx J dy ~p2 i 2 * 2 ™* + - xm v) [(“£ + bd ) + ( bd - ac ) 75 ] 

0 0 


On the other hand the contribution of the vacuum polarization vanishes. Therefore, the full contribution to the MDM 
yields 

A 2 HDM (<h 1 )mdM = A 2 H±IL (?’ 1 )mDM + A 2Lli/± 1 )mDM ( 22 ) 

for the 2HDM type I, the values of a, b , c and d are 


a = c = 


tan/3 mvi 


b=d = 


tan/3 


rtu 


we shall use the numerical value _ 

G f = = 1 1663787(6) x 10- 5 GeK- 2 

as for the 2HDM type II, the values of a,b,c and d are given by 


a = c = 


2 tynj 

tan 0 


■m Vl 


b,d = 2* y/G ftyii tan /3 


6 Results and analysis 

Our analysis will be based on constraints on charged Higgs masses and the tan /3 parameter. For either model type I or 
II the experimental constraints on the possible values in the (m.H± > tan 3) parameter space comes from processes such as 
B u —> rv T , D s —> r^ T , B —» Dtu t , K —> pv^ and BR (B —y X.7)[I5]. 

Based on the phenomenological constraints on the 2HDM type I, we take values of tan/3 between (2 — 90) and values 
of the charged Higgs mass of mn± = (100 — 300 — 500 — 700 — 900) GeV |17|. On the other hand, for the 2IIDM type 
U, we have different allowed intervals of tan/3 for different values of the charged Higgs mass: for m//± = 300 GeV 
the values of tan/3 lie within the interval (4 — 40), for m//± = 500 GeV the value of tan/3 is between (2 — 69) and for 
m H ± (700 — 900) GeV the values of tan/3 is between (1 — 70) [17] . 

We shall make contourplots of the neutrino mass versus MDM of the neutrino for different values of the charged Higgs 
mass sweeping all allowed values of tan/3 for each mass. As for the neutrino masses, we shall plot up to an order of 
magnitude higher than the upper bound of the SM. 

• Electron neutrino case 

Taking into account the upper experimental bound in the SM for the electron neutrino mass rn Ue , we shall plot within 
the interval 1 x 10 ~ 8 MeV < m Ue < 1 x 10 -5 MeF. If we use the above interval in a model with two Higgs doublets, 
for different values of Higgs mass and tan/3, we shall obtain exclusion regions by taking as reference the experimental 
thresholds for the MDM of the electron neutrino. In this way, it is possible to obtain upper bomids on the neutrino mass 
in this scenario. 

In Fig. 0 we plot the electron neutrino mass versus MDM for charged Higgs masses of 100, 300, 500, 700,900 GeV for 
the 2HDM type I (left-hand side) and for masses of 300,500,700,900 GeV for the 2HDM type II (right-hand side). The 
horizontal lines correspond to the experimental upper limits for MDM coming from TEXONO 2007 (Taiwan Experiment 
On NeutriNO) [18] which is p Ve < 7.4 x 10 -1 Vb at 90% C.L. , and GEMMA 2013. (Germanium Experiment for 
measurement of Magnetic Moment of Antineutrino)[19) which is p„ e < 2.9 x 10 -11 /iB at 90% C.L.. We observe that the 
maximum values of MDM that can be reached are 1.635 x 10 ~ 17 Pb f° r a value of the charged Higgs mass of niff± = lOOGeV' 
and tan/3 = 2 in the case of the 2HDM type I with an electron neutrino mass of 1 x 10 ~ 5 MeV and 1.715 x 10 ~ 17 pb for 
a value of the charged Higgs mass of m//± = 500 GeV and tan/3 = 69 in the case of the 2HDM type II with an electron 
neutrino mass of 1 x 10 ~ b MeV . These values are far from the experimental threshold and provide no bounds on the 
neutrino mass. 
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Figure 6: The graphics show the values of the magnetic dipole moment (MDM) as a function of the electron neutrino mass 
between (l x 10 -8 -lx 10 -5 ) MeV and for values of the charged Higgs mass of (100 - 300 - 500 - 700 - 900) GeV for 
the 2HDM type I (left-hand side) and (300 — 500 — 700 — 900) GeV for the 2HDM type II (right-hand side). 


• Muon neutrino case 


We shall plot within the interval 2 x 10 ~ 6 MeV < rn Utl < 4 x 10 ~ 1 MeV. Using such an interval in the 2HDM, 
for different values of Higgs mass and tan /3, we can obtain exclusion regions by taking as reference the experimental 
thresholds for the MDM of the muon neutrino. In that way we obtain upper limits for the muon neutrino mass in this 
scenario. 


1E . 9 MDM Ipg > 
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Figure 7: Values of the magnetic dipole moment (MDM) as a function of the muon neutrino mass between 
(1 x 10 -5 - 4 x 10 -1 ) MeV and for values of the charged Higgs mass of (100 - 300 - 500 - 700 - 900) GeV for the 2HDM 
type I and (300 — 500 — 700 — 900) GeV for the 2HDM type II. 


In Fig. 0 we plot the muon neutrino mass versus MDM for the same charged Higgs masses as before for the 2HDM 
type I (left-hand side) and type II (right-hand side). The horizontal lines correspond to the experimental limits for MDM 
coming from LSND 2001 (Liquid Scintillating Neutrino Detector) [20] which is /i„ M < 6.8 x 10~ 10 ^b 90% C.L., and 

BOREXino 2008 (Boron solar neutrino experiment) pH] which is p Ult < 1.9 x 10” 1 ( Vb at 90% C.L.. 

From Fig. |7] we can see that it is necessary to make a more detailed analysis for certain values of the Higgs mass 
and the respective values of tan/3 for each model, in order to find upper limits for the neutrino mass as a function of 
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Figure 8: Values of the MDM as a function of the muon neutrino mass between (lx 10 2 — 4 x 10 *) MeV and for 
values of the charged Higgs mass of 100 GeV for the 2HDM type I. 


tan/3 

™>H± 

GeV) 

Experiment 

300 

500 

700 

900 

2 

- 

- 

- 

- 

LSND 2001 

— 

3.31 x 10" 1 

— 

— 

BOHEXIiNO 2008 

4 

1.85 x 10" 1 

— 

— 

— 

LSND 2001 

6.47 x 10 z 

— 

— 

— 

BOREXINO 2008 
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1.16 x 10" A 

2.82 x 10' 1 

— 

- 

LSND 2001 

5.3 x 10 z 

8.71 x 10 z 

2.18 x 10" 1 

2.25 x 10~ ‘ 

BOREXINO 2008 

15 

3.95 x 10"* 

6.16 x 10~ z 

8.64 x 10~ z 

1.16 x 10” 1 

LSND 2001 

2.07 x 10~ z 

3.18 x 10 * 

4.26 x 10" z 

5.27 x 10~ z 

BOREXINO 2008 

30 

2.2 x 10 z 

3.36 x 10--= 

4.51 x 10"-* 

5.57 x 10"-* 

LSND 2001 

1.15 x 10 z 

1.78 x 10 z 

2.39 x 10 z 

2.89 x 10 z 

BOREXINO 2008 

40 

1.72 x 10"* 

— 

- 

- 

LSND 2001 

8.95 x 10 J 

- 

— 

— 

BOREXINO 2008 

69 

— 

1.67 x 10“ z 

— 

- 

LSND 2001 

— 

8.7 x 10 J 

— 

— 

BOREXINO 2008 

70 

— 

— 

2.19 x 10~ z 

2.69 x 10" z 

LSND 2001 

— 

— 

1.15 x 10~ z 

1.43 x 10" z 

BOREXINO 2008 


Table 3: This table shows upper bounds for the muon neutrino mass (MeV) as a function the free parameters tan 0 and 
rtiH± in the 2HDM type II, taken from figure [9] The empty cases correspond to excluded regions of the model. 
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Figure 9: Values of the MDM as a function of the muon neutrino mass between (l x 10 3 — 4 x 10 *) MeV and for 
values of the charged Higgs mass of (300 - 500 - 700 - 900) GeV to 2HDM type II 
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charged Higgs masses, tail /3 and the current experimental limits. This more detailed analysis is shown in Fig. [H] for the 
2HDM type I and in Fig. [T7| for the 2HDM type II. We observe that in the 2HDM type I the strongest bound for the 
muon neutrino mass that we obtain is given by 2.583 x 10~ 1 MeV for M^+ = 100 GeV and tan/3 = 3 obtained from the 
bound of MDM from BOREXino. Nevertheless, for the interval of neutrino mass plotted we do not obtain bounds on the 
neutrino mass from the bound of MDM coming from LSND, neither for other masses of the charged Higgs. 

As for the 2HDM type II, we observe significant differences in the patterns of the bounds. For instance, in the 2HDM 
type II the hierarchy of the bounds are in opposite order as a function of tan/3 with respect to the 2HDM type I. This 
happens because in the 2HDM type I the couplings are proportional only to tan/3, while for the 2HDM type II the 
couplings are proportional to either tan /3 or cot /3. Of course, when the charged Higgs mass increases, the bounds on the 
muon neutrino mass become less restrictive since the contribution of the new physics tends to decouple as the Higgs mass 
increases. 

Considering the current experimental limits, the strongest upper limit that could take the muon neutrino mass in the 
2HDM type II is 8.953 x !0~ 3 MeV, obtained from the set of parameters m//± = 300(jeV and tan/3 = 40. The numerical 
values of the upper bounds for the muon neutrino mass within the 2HDM type H, are shown in table 0 for different 
allowed values of ra//± and tan/3 parameters. 
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Figure 10: Values of the MDM as a function of the tau neutrino mass between (lx 10 6 — 3 x 10 1 ) MeV and for values 
of the charged Higgs mass of (100 - 300 - 500 - 700 - 900) GeV for the 2HDM type I and (300 - 500 - 700 - 900) GeV 
for the 2HDM type II. 


• Tau neutrino case 

We shall plot within the interval 2 x 10 -6 < m Ur < 20A leV, and obtain our bounds from the experimental limits on 
the MDM. In Fig. [TTB we plot the tau neutrino mass versus MDM for the same charged Higgs masses as before for the 
2HDM type I (left-hand side) and type II (right-hand side). The horizontal lines correspond to the experimental limits 
for MDM coming from DONUT 2001(Direct Observation of the NU Tau)[22] which is p l , T < 3.9 x 10 -7 /zb at 90% C.L., 
and BOREXino 2008[2T] which is /i„ M < 1.5 x 10 ~ 10 pb at 90% C.L. 

Figure [TD] shows that we require a more detailed analysis of the upper bounds of the tau neutrino masses in terms of 
the free parameters. Such an analysis is carried out in Fig. [IT] for the 2HDM type I and in Fig. [T2]for the 2HDM type II. 
Once again we have significant differences in the patterns of the bounds for the models of type I and of type II because of 
the different behavior of the couplings with respect to the tan/3 parameter. Further, the upper limits of the tau neutrino 
mass is weakened as the charged Higgs mass increases owing to the decoupling behavior of the diagrams with respect to 
the Higgs mass. 

Considering the current experimental limits, the strongest limit obtained for the tau neutrino mass in the 2HDM type 
I is 9.098 x 10 _2 A/eV, and occurs for the values m H ± = 100 GeV and tan/3 = 3. As for the 2HDM type H the strongest 
upper limit on the tau neutrino mass is 2.926 x 10 -3 A/eV obtained with the set of parameters m//± = 300 GeV and 
tan /3 = 40. Finally, the numerical values of the upper bounds for the muon neutrino mass within the 2HDM type I and 
II, are shown in tables0and|]3respectively, for different allowed values of m//± and tan/3 parameters. 
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Figure 11: Values of the MDM as a function of the tau neutrino mass between (lx 10 2 — 3 x 10 1 ) MeV and for values 
of the charged Higgs mass of (100 - 300 - 500 - 700 - 900) GeV for the 2HDM type I. 
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Figure 12: Values of the MDM as a function of the tau neutrino mass between (lx 10 3 — 3 x 10 1 ) MeV and for values 
of the charged Higgs mass of (300 - 500 - 700 - 900) GeV for the 2HDM type II 
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tan/3 

m H ± (GeV) 

Experiment 

100 

300 

500 

700 

900 

2 

- 

10.95 

17.42 

23.95 

29.78 

DONUT 2001 

BOREXINO 2008 

— 

1.55 

2.71 x 10" 1 

3.95 x 10“ 1 

5.3 x 10“‘ 

3 

6.28 

— 

— 

— 

— 

DONUT 2001 

BOREXINO 2008 

8.86 x 10" z 

— 

— 

— 

— 

7 

16.67 

- 

- 

- 

- 

DONUT 2001 

BOREXINO 2008 

2.6 x 10" 1 

7.08 x 10 1 

1.19 

1.74 

2.19 

20 

- 

— 

— 

— 

- 

DONUT 2001 

BOREXINO 2008 

7.74 x nr 1 

2.31 

3.75 

5.29 

6.59 

40 

— 

— 

— 

— 

— 

DONUT 2001 

BOREXINO 2008 

1.72 

4.84 

7.69 

10.4 

13.07 

60 

— 

— 

— 

— 

— 

DONUT 2001 

BOREXINO 2008 

2.67 

7.19 

11.81 

15.68 

19.35 

90 

— 

— 

- 

- 

- 

DONUT 2001 

BOREXINO 2008 

4.05 

10.95 

17.42 

23.64 

29.4 


Table 4: This table shows the upper bounds for the tau neutrino mass ( MeV ) for several allowed values of the free 
parameters tan/3 and rrin± in the 2HDM type I, taken from figure [TT] The empty cases correspond to excluded regions 
of the model. 


tan/3 

m. H ± 

GeV) 

Experiment 

300 

500 

700 

900 

1 

— 

— 

11.87 

15.65 

DONUT 2001 

— 

— 

1.73 x 10" 1 

2.27 x 10-* 

BOREXINO 2008 

2 

- 

4.44 

- 

- 

DONUT 2001 

— 

6.24 x 10 “ 

— 

— 

BOREXINO 2008 

4 

1.24 

— 

— 

— 

DONUT 2001 

2.27 x 10 z 

— 

— 

— 

BOREXINO 2008 

5 

9.52 x HT 1 

1.58 

2.27 

2.97 

DONUT 2001 

1.82 x 10 z 

2.77 x 10 “ 

3.66 x 10~ z 

4.52 x 10~ z 

BOREXINO 2008 

15 

2.48 x 10' 1 

4.58 x 10~ A 

6.49 x 10- 1 

8.63 x 10- 1 

DONUT 2001 

7.13 x 10 J 

1.13 x 10 “ 

1.43 x 10~ z 

1.73 x 10~ z 

BOREXINO 2008 

30 

1.14 x 10' 1 

1.93 x 10 1 

2.89 x 10" 1 

3.82 x 10“ l 

DONUT 2001 

3.7 x 10 J 

5.82 x 10 J 

8.04 x 10 J 

9.62 x 10 ° 

BOREXINO 2008 

40 

8.58 x 10 z 

— 

— 

— 

DONUT 2001 

2.87 x 10 0 

— 

— 

— 

BOREXINO 2008 

69 

— 

8.14 x 10’^ 

— 

— 

DONUT 2001 

— 

2.71 x 10 J 

— 

— 

BOREXINO 2008 

70 

— 

— 

1.12 x 10~‘ 

1.39 x 10“ 1 

DONUT 2001 

— 

— 

3.57 x 10-° 

4.48 x 10-° 

BOREXINO 2008 


Table 5: Upper bounds for the tau neutrino mass (MeV) for several allowed values of the free parameters tan/3 and m//± 
in the 2HDM type II, taken from figure [T2] The empty cases correspond to excluded regions of the model. 
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7 Conclusions 

The neutrino magnetic moment provides a tool for exploration of physics beyond the Standard Model. The magnitude 
of the magnetic moment is highly sensitive to the neutrino mass, but also depends on the mass of the associated charged 
lepton inserted into the loops. 

Further, the value of the magnetic moments of the neutrinos could be modified with physics beyond the Standard 
Model, hi particular for the Two Higgs Doublet Model (2HDM) we evaluated the contributions coming from the insertion 
of the charged Higgs boson into the loops. Our results show that for the 2HDM of type I and of type II, the total 
contribution is far from the threshold of experimental detection in the case of electron neutrinos (owing to the supression 
coming from the electron mass into the loops), obtaining a maximum contribution about six orders of magnitude below 
the present experimental limits, hi the case of muon neutrinos the total contribution produces weak bounds for the mass 
of the neutrinos for model type I, and stronger bounds for the case of model type II. Finally, such bounds are much 
stronger for tau neutrinos (because of the enhancement of the tau mass into the loops) for either type of 2HDM, but 
restrictions are much stronger for the model type II. 

hi general since the bounds are highly sensitive to the value of the tan ft parameter, the limits obtained are significantly 
different for the model type I with respect to the model type II because of the different dependence on the Yukawa 
couplings of each model with tan/?. Of course, the limits are weakened as the mass of the charged Higgs increases 
since the contribution of new physics tends to decouple as the Higss mass grows. Further, since the contribution of 
diagrams involving the associated charged lepton increases with the mass of the charged lepton, the strongest bounds are 
obtained for the neutrino associated with the heaviest charged lepton (tau neutrino) while basically no bounds near the 
experimental threshold are obtained for the electron neutrino. 
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A Explicit expressions for the EFF 

In this Appendix we present some details of the process of calculating the EFF’s. For the case of two higgs bosons and 
one lepton (2H ± 1L), the general form of the contribution can be written as 

df'k \ (2 k a + p? + pf) (aP L + bP R ) (/l + m,) (cP L + dP R ) 


, ,, f d*k f ( 2k° 


+ Pi) 2 — m //±] [(k + P 2) 2 — m %± (fc 2 - mf) 
expanding the numerator, denoting A = (aP L + bP R ) and D = (cP L + dP R ) we find 

(2 k a +p°+pf)A(lk + m,)B = (2 k°kg + p^k B + p°k g ) A-, a B + m, ( 2k° + pf + pf) AB 
and for the denominator we use the dimensional regularization method 

1 


(fc+P2> - m 2 H± \(k 2 - mf) (A; + pi) - m 2 H ± 


ai 


a 2 


03 


then 


a } a 2 a 3 


r(l + l + l) 

r (i) r (i) r (i) 

1 x 

= r(3) J dx j dy 


hh 

0 0 


' [a 3 (x 0 -x)+a 2 (x-y) + a\ (y - z 3 )] 

1 


0 0 [((fc + pi) 2 - m 2 H± )(!-«) + (fc 2 - mf) {x-y)+ [(fc + p 2 ) 2 - m 


2 

H± 


~\ n 3 
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where xo = 1 and X 3 = 0. Thus, the denominator can be written as 

[(k + ptf - m 2 H± {l-x) +(k 2 -ni})(x-y)+[(k+ p 2 ) 2 -m 2 H± y 
= k 2 + 2k- (pi (1 - x) + P 2 y) + (m%± - ml- m 2 ) x + (m? + m 2 - m 2 H ±) y + m 2 - m 2 H ± 

'-*-' ^--—v-' 

b a 2 

= k 2 + 2k - b + a 2 
= k 2 + P 2 

where we use the transformation k —» k — b in the last equation and 

P 2 = a 2 -b 2 = y 2 ml - 2xym 2 u + (m 2 - m} + m 2 H± ) y + x 2 m 2 u + (mf - ml - m 2 H± ) x + m 2 H ± 
in consequence, adding the corresponding terms to the integral with terms in the numerator 1 , and k^k^ we obtain 

cPk ( (2k a ks + pf ka + p? kp) A^B + mi (2k° + p 2 + pf) AB 


^2H±1L 


™ - -m{ 


(k + pi ) 2 - m 2 H± j [(k + P 2) 2 - m 2 H± J (k 2 - mf) 
2 b Q b 0 


= —e 


& J */% v- + ( 


26“ _ p2a _ plo 
p2 p2 p2 


0 0 

expanding the terms bP and employing the Dirac equation 


miAB 


(7 ■Pi-m)u(p 1 ) = 0=^u(p!) = mu(p!) 
w(p 2 )( 7 ’P 2 ~m) = 0 =»u(p 2 )^ 2 = mu(p 2 ) 


then 


A?*±u .(®,0 


= —e 


= —e 


16tt 2 


16tt 2 


1 X 

/' lx / rfy [( 


26 “ frS , c$, A2 , P2a& a , Pl a b ff \ t .. B D , f26“ P 2 a Pla 1 „ „ 

-pT-+9 a hp2+-p2- + ^pT-) A r D +{-pT- ^2~ -p2\ m > AD 


0 0 

hi 

0 0 


dy 


1 

p2 


-Pia (1 - x) 2 - p 2a y (1 - x) - p la (1 -x)y- p 2a y : 
+ (p 2 a (1 - x) + p 2a y) £ + (Pla (1 — ar) + Pla2/) ? 


m„ ((ac + bd) + (bd — ac ) 75 ) 


I . v 11 

+m I Pla (1 - X) + p 2a y - P 2 c- - Pla- 


[ac + bd) + (bd - ac) 75 ] J + - 7 ° [( 6 c + ad) - (be - ad) 75 ] 


-(£)] 


and using the Gordon relation 


1 


ft (P 2 a) 7 aU (pia) = (p 2 a) [(P2 + Pl) Q + w (Pla) 

=> ft(p2o) (P2 + Pl) a U (pi Q ) = 77 (p 2a ) 2m v ^ Q U (pi Q ) - ft(p2a) UT^q^U (pla) 

finally the contribution for the EFF’s with two charged Higgses and one lepton can be represented by 
A? 


L 2H±1L (^*0 “ 
1 x 

—ie 


J dx f dy ^2 ( ( ( m/m " ” 2 m v) 7a + X (6ml - 2m„mi) 7 Q - 4x 2 m 2 7 Q + ^ 


0 0 


mi 1 io Qtl qn 


(m t - 3m„) ixo otl q ll + 2x 2 im l/ o Qtl q ll ) [(ac + bd) + ( 6 d - ac) 75 ] + ^ 7 “ [( 6 c + ad) - (be- ad) 75 ] hi ^^ 2 ^ 
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And for the case of two lepton and a charged Higgs ( 2L1H ± ) the general form of the contribution can be written as 

f d 4 k ( (“-Pi. + bPn) (* + »*!+ m ‘) 7“ {t + f 2 + mi) ( cP L + dP R ) | 

;WM) = - e y ^ / 


expanding the numerator and using the same change as for the vertex 2H ± 1L we have 

A (# + h + m i) 7° + + m 0 B = 

+ pi^kp + k li p 20 + P\tf> 2 p) A'f'f'fB + mi (k + A^°B + m t (k + p 2 ) 0 A'f'fB + mfAy a B 
and for the denominator we use the dimensional regularization method 
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(k + p 2 ) - mf (k 2 - m 2 H ±) (A; + pi) - mf 


“1 


«3 


then 




[((fc+pi ) 2 - mf) (l-x) + (k 2 — mjj±) (x-y) + ((k + p 2 f -mf) y 
where xq = 1 and £3 = 0. Hence, the denominator can be written as 


[(A: + pi ) 2 - mf (1 - x) + (A : 2 - m 2 H ±) (x - y) + [(A; + P 2) 2 - 


mf 


= A : 2 + 2k • (pi (1 - x) + p 2 y) + (mf - mf - m 2 H ±) x + (m 2 H ± + mf - mf) y + mf - mf 

6 a* 

= A : 2 + 2k • 6 + a 2 
= A ; 2 + P 2 

where we use the transformation k —> k — b in the last equation and 

P 2 = a 2 -b 2 = y 2 m 2 - (m 2 H ± - mf - mf) p - 2xym 2 „ + x 2 mf - (mf + mf - m 2 H± ) x + mf 
therefore, adding the corresponding terms to the integral with terms 1,A; / * and k^k^ in the nmnerator, we obtain 


= — e 
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expanding the terms fe M and employing the Dirac equation we obtain 


^2L1H± (fl»0 



0 


0 0 


0 


+ (£ - S) m ‘ A ^ B - 


1 X 



0 0 


0 0 


+ ((y - x + xy) (P 2 a + pio) - p2 a y 2 +Pia(x-y — x 2 )) m„ ((ac + fed) + (fed - ac) 75 ) 
m/m*, 7 ° ((fee + ad) - (fee - ad) 75 ) + mi (p 2 a (p - 1 ) - zpia) ((ac + fed) + (fed - ac) 75 ) 


((fee + ad) - (fee - ad) 75 ) J ~ -^7° ((^ + ad) - (fee - ad) 75 ) In ^ 


((fee + ad) - (fee-ad) 75 ) 


and using the Gordon relation like in the previous case, the contribution to the EFF’s with two leptons and one charged 
Higgs can be represented as 


^2L1H± (tf*0 — 


1 X 



d) - (fee - ad) 75 ] ln-^- 
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